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We introduce and study an extended "t-U-J" two-orbital model for the pnictides that includes 
Heisenberg terms deduced from the strong coupling expansion. Including these J terms explicitly 
allows us to enhance the strength of the (tt, 0)-(0, tt) spin order which favors the presence of tightly 
bound pairing states even in the small clusters that are here exactly diagonalized. The Aig and 
B2g pairing symmetries are found to compete in the realistic spin-ordered and metallic regime. 
The dynamical pairing susceptibility additionally unveils low-lying Big states, suggesting that small 
changes in parameters may render any of the three channels stable. 



Introduction. One of the main puzzles in iron-based su- 
perconductors [l| arises from the conflicting experimen- 
tal results on the presence of nodes in the superconduct- 
ing state. Surface-sensitive angle-resolved photoemission 
(ARPES) studies indicate that full nearly momentum- 
independent gaps open on all Fermi surface (PS) pock- 
ets. However, some bulk experimentsgive results com- 
patible with nodal superconductivity |3j. On the theory 
side, calculations where many different pairing states are 
allowed to compete, as opposed to studying a few iso- 
lated states, are difficult for multiorbital Hubbard mod- 
els. Within magnetic mechanisms for superconductivity, 
two approaches have addressed this issue: (i) random 
phase approximation (RPA) studies suggested that sev- 
eral pairing channels are in competition [j| , in agreement 
with (ii) two-orbital model Lanczos studies [5|,@ based on 
the quantum numbers of the state with two more elec- 
trons than half-filling. However, RPA relies on a par- 
ticular subset of diagrams, while studying the quantum 
numbers of clusters did work before for the cuprates , 
but finding true Cooper pair formation is difficult. 

In this Letter, a simple generalization of Hubbard mod- 
els for pnictides is presented that increases the strength 
of the (tt, 0)-(0, tt) spin order in the undoped limit, cre- 
ating tightly bound-states upon electronic doping that 
can be studied with Lanczos methods on the small clus- 
ters currently accessible with state-of-the-art computers. 
These extra terms can be justified by noting that the 
undoped state, which combines itinerant electrons with 
a robust Neel temperature and a lattice distortion Q, 
is itself rather exotic and suggests a more stable mag- 
netic order than a description based exclusively on on- 
site Coulomb repulsion would support. The t-J model 
for cuprates provides further guidance: here J, when 
considered as independent of t, can be increased to suf- 



ficiently large values that d-wave pairing tendencies are 
amplified and tightly bound-states are formed, while in 
the Hubbard model the d-wave pairing signal is weak 0] . 
In pnictides, the strong coupling i- J1-J2 model has been 
studied before [l^. In the alternative "t-[/-J" route [llj 
to be followed here, Heisenberg "J" terms will be added 
to the original Hubbard model to enhance spin order and 
pairing tendencies, but without projecting out doubly oc- 
cupied sites and charge fluctuations. 

Model and method. The two-orbital model 0, H, 
based on the dxz {x) and dyz (y) Fe orbitals is studied 
here. Keeping only these two orbitals is reasonable since 
X and y provide the largest contribution to the pnictides' 
band structure PS • In addition, the studies described 
below are computationally demanding and they simply 
cannot be carried out with more orbitals. Thus, a bal- 
ance must be reached between the more ideal five-orbital 
models and the feasibility of the actual calculations. The 
model includes a hopping term with amplitudes that fit 
band calculations [12| (energy scale \ti\), a Hubbard term 
with on-site intraorbital repulsion U, a Hund coupling 
Jh, an interorbital repulsion fixed as U' = J7 — 2Jh, 
and a pair-hopping term with coupling J'=Ju [l3l- The 
model was used in several previous studies [61 . The nov- 
elty are the extra Heisenberg terms that will be added, 
with associated exchange couplings for nearest-neighbor 
(NN) and next-nearest-neighbor (NNN), Jnn and Jnnn, 
as discussed below. 

The two-orbital Hamiltonian is exactly investigated us- 
ing the Lanczos algorithm (including dynamical informa- 
tion) on a small tilted VSxVS cluster 0, 0,11]. This re- 
quires substantial computational resources: to determine 
the undoped-limit ground state of the 8-sites cluster even 
exploiting the Hamiltonian symmetries still requires a ba- 
sis with ~2-20 M states (equivalent to a 16-sites cluster 



2 



one-band Hubbard model) depending on the subspaee ex- 
plored. Lanczos runs had to be performed for all allowed 
momenta k, quantum numbers under rotations and re- 
flections (i.e. irreducible representations Aig, A2g, Big, 
i?2g, and E of the Z?4h group @, and z-axis total 

spin projections. In addition, binding energies require 
calculations for a number of electrons A*" equal to 16, 17, 
and 18, varying U, Jh, Jnn, and Jnnn in a fine grid. 
The full effort amounted to ^8,000 diagonalizations of 
the cluster, supplemented by dynamical calculations, us- 
ing a Penguin 128GB Altus 3600 computer. 

Results for the original two-orbital model. The relative 
symmetry between the undoped (iV=16) ground state 
(GS) and the N=18 GS has been studied varying U/\ti\ 
and Jb_/U. The results at Jnn=>/nnn=0 are shown in 
Fig. [TJa). The undoped GS has momentum k = (0,0) 
and it transforms according to the Aig representation of 
the D^h group, for all the investigated values of Jh and 
U. The 7V=18 GS also has k = (0,0), but its irreducible 
representation varies in different regions of the phase di- 
agram. In agreement with previous results the iV=18 
GS is a spin triplet for U< 6 ^^'^ ^ broad range 

oi Jk/U. a spin-singlet state with symmetry B2g dom- 
inates the small Ju/U (roughly <0.15) region for the 
studied values of U /\ti\. This B2g state arises from the 
orbital portion of the pairing operator since the two extra 
electrons added to the undoped GS are located in differ- 
ent orbitals At C/> 7 and intermediate to large 
Ju/U regimes, the singlet-state symmetry becomes Aig. 
The binding energy Eb=E{18) + E{16) - 2£'(17), where 
E{N) is the GS energy of N electrons, was also studied. 
Binding, i.e. Eb < 0, is observed but only at large J7's 
where the undoped GS is an insulator. 
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FIG. 1: (color online) Relative symmetry between the A*'=16 
(undoped) and A'^=18 GS's, varying U and Jh/U. Circles 
denote triplet states, squares _B2g-symmetric singlets, and di- 
amonds Aig-symmetric singlets, (a) Results for couplings 
Jnn= Jnnn = 0. Open triangles indicate binding, (b) Results 
for the lowest value of (Jnn, Jnnn) where binding appears. 

Binding stabilization. As discussed in the Introduc- 
tion, here the spin background with wavevectors (tt, 0)- 



(0, tt) will be magnified via the addition of extra Heisen- 
berg terms with the expectation that carrier attraction 
will become stronger, leading to Eb < pairing. To 
find the precise form of the Heisenberg terms, a strong 
coupling (large U and Jh) expansion of the undoped 
two-orbital Hamiltonian was carried out. NN and NNN 
Heisenberg interactions Jnn and Jnnn were obtained. 
This strong coupling expansion for multiorbital models 
is subtle and it has to be performed in such a way that 
the results are independent of the basis chosen for the or- 
bitals. 
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sites with J,-, 



these additional interactions 
where (ij) indicates NN and NNN 
Jnn and Jnnn, {q:,/3} label the orbitals 



(x or y), and Sia is the total spin at site i and orbital 



and Jnnn= 



2(t -\-t ) 

a. The couplings are Jnn= 3([/^j^-) 
where ti are the hoppings ti=—l, t2~1.3 
^ i\ti \ units) 



3{U+Jh) ' 

i3=i4=-0.85 
Jnn and Jnnn are 
^-^$^=0.93 
antiferro- 



in the usual notation 
the same for both orbitals. The ratio -r- 

-^NNN 

is kept fixed since Jnn and Jnnn are both 
magnetic, and a ratio ~1 introduces frustration favoring 
(tt, 0)-(0, tt) spin order, as confirmed by calculating the 
spin structure factor and varying Jnn/C^ [13 ■ Below, 
Jnn will be considered a free parameter independent of 
U and Jh (with Jnn/ Jnnn=0.93), to further enhance 
such a spin order. 

By adding the extra Heisenberg terms to the two- 
orbital model, the desired goal is obtained since with 
increasing Jnn eventually Eb becomes negative for all 
the studied ([/, Jh) couplings. The spin-triplet region 
virtually disappears (Fig. [IJb)) and it is mainly re- 
placed by the i?2g state which itself becomes confined 
to U < 4 1^1 1 (squares) due to the expansion of the Aig 
region. The symmetries shown in Fig.[TJb) were obtained 
with the smallest superexchange values ( Jnn' "^nnn) that 
produce binding of two electrons at each ([/, Jh) point. 
In Fig. [2l^a), where the binding energy vs. J]<i^/U is 
shown at several J7's and at a fixed (realistic) Jh/J7=0.2, 
some examples of (Jnn, '^nnn) found. Increasing 

Jnn eventually induces binding for all J7's. The value of 
Jnn/U for which binding occurs decreases as U increases. 

Studying Eb and the relative symmetry between the 
A''=16 and 18 GS's, phase diagrams in the (J7, Jnn/C/) 
plane were constructed. In Fig. [2ljb), typical results for 
Jh/?7=0.2 are shown [l3|. The bound state has Aig sym- 
metry in most of the binding region, but a i?2g symmet- 
ric state also prevails at small U values (~ 2 \ti\). Both 
symmetries appear inside the proper magnetic/metallic 
region of the undoped limit, according to mean-field cal- 
culations [i^l extended to incorporate Jnn- In [13 it is 
shown that the results in Figs. [2l[a,b) arc qualitatively 
the same varying Jnn/Jnnn in the range [0.5,1.5]. 

Overlaps. Consider now the pairing operators that 
produce the electronic bound states. The over- 
lap (^'(iV=18)|Aj^ J*(7V=16)) was calculated, where 
|5'(A^)) is the GS in the subspaee of electrons and 



FIG. 2: (color online) (a) -Be/ltil vs. Jnn/I/ for different 
values of (7/|ti| and Jh/J7=0.2. (b) Phase diagram show- 
ing "Binding" and "No Binding" regions and the symmetry 
of the two-electron bound state varying t//|ti| and Jnn/U, 
for Jh/U=0.2. The shaded area is where the antiferromag- 
netic/metallic state is stabilized in the mean- field approxima- 
tion for the undoped limit. The dot-dashed line is for Fig. O 



^k,i=Sa/3/W('^')"f4,a,T4.^.i' ^'^^ '^k.a,. Creating 
an electron with spin z-axis projection a, at orbital 

a = x^y, and with momentum k. The structure factor 
/(k) arises from the spatial location of the electrons form- 
ing the pair and Gi are the Pauli matrices {i = 1, 2, 3) 
or the 2x2 identity matrix cto (* = 0) (note that cri and 
CT2 imply an interorbital pairing). Overlaps for all the 
symmetries in and with NN and NNN locations for 
the electronic pairs, were evaluated. 

For all operators respecting the relative symmetry be- 
tween the doped and undoped states, finite overlaps 
were found, although of different values. As a trend, 
as the binding grows, pairing involving NNN opera- 
tors prevail over the NN ones. For example, in the 
Aig region in Fig. [^b) there arc four pairing opera- 
tors with finite overlap (shown in Fig. for [/=3|ti| 
and Jh/?7=0.2) characterized by /(k)cri equal to: (i) 
(cos fcj: 4- cos fcj,)cro (full circles); (ii) (cos /cj. cos fcj,)cro (full 
squares); (iii) (sin fc^; sin fcy)cri (full diamonds); and (iv) 
(cosfca; — cosfcj,)cr3 (full triangles). Close to the bound- 
ary with the phase where the binding is weak [E-qp^- 
0.05 |ii|), operators (i) and (ii) present the largest, and 
almost equal, overlaps. With increasing binding the (i) 
overlap decreases while (ii) becomes stronger. The over- 
laps for operators (iii) and (iv) are clearly smaller. 

Note that (ii) is the simplest expression of a nodeless 
s± pairing operator [2l|. Our results indicate that this 
type of pairing dominates only when the binding energy 
is large, which occurs at very large U or Jnn- At interme- 
diate values of couplings, a symmetric linear combination 
of (i) and (ii) with almost equal weights is optimal, and 
it leads to a "quasi-nodal" s± pairing state (Fig. [3Jb)). 
From this perspective, the most "natural" Aig pairing 



FIG. 3: (color online) (a) Overlap (*(iV=18)|Aj, .|*(Ar=16)) 
vs. Jt^j-i/U for the indicated pairing operators, at [7=3 \ti\ 
and Jh/U—0.2. (b) Superconducting gap at the FS: internal 
hole pocket (continuous line), external hole pocket (dashed 
line). The dot-dashed and double dot-dashed lines are for 
the two electron pockets which intersect at the Brillouin zone 
boundary ( $=7r/4) of the folded zone. The Aig+ symmetric 
linear combination of Aig operators (i) and (ii) is used, with 
equal weight. The angle $ is measured from the positive x- 
axis to the positive y-axis. (c) Same as (b) but for the -B2g+ 
symmetric combination of the B2g operators (v) and (vi). (d) 
Dynamic pairing susceptibility for the pairing operators indi- 
cated (see text), at (7=3 Jn/U^O.2, and Jnn/(7=0.095. 
The vertical line indicates -B(18) — £(16). 

operator arises from a linear combination of (i) and (ii), 
as opposed to just (ii) as in s± scenarios. The gaps in 
Fig-ETb) were calculated from mean-field approximations 
as in and 2^ , and choosing a pairing strength Vq such 
that the gap order-of-magnitude in meV's agrees with ex- 
periments [16| . Note that the linear combination Aig+ 
for the hole pockets closely reproduces (full and dashed 
lines) the ARPES results in the superconducting state, 
with both gaps only weakly k-dependcnt, and with the 
interior (exterior) pocket gap ^ 12 (6) meV. The electron 
pockets, on the other hand, present strongly k-dependent 
gaps, and a quasi-node is found along the x- (y-) axes for 
the pocket at X (Y). In the folded zone, this implies that 
the quasi node is on the outer pocket, in a gree ment with 
angle-resolved specific heat measurements (23[ . 

Note that the presence of a dxy "patch" on the elec- 
tron pockets has been discussed before by many groups 
as possibly responsible for gap nodes (or minima) on the 
electron pockets. The present results show that such a 
minimum (or nodes) can arise without such an xy-patch, 
which is important to assess the impact of the various 
orbitals. The one-particle spectral function ^(k, w) was 



also calculated [18| . Features on the scale of the magnetic 



or superconducting gaps cannot be resolved within the 
few momenta available, but the higher energy features 
at intermediate couplings arc similar to non-interacting 
bands [20| . in agreement with ARPES experiments and 
with local density plus dynamical mean-field theory, cal- 



4 



culations |24l |. 

As mentioned before, in physically relevant portions of 
the phase diagram [2^ the pairing symmetry com- 
petes with Aig. Three pairing operators with finite 
overlaps were fomid in this region: (v) (cosfca;-l-cosfc.y)(Ti, 

(vi) (cos kx cos ky)a I, and (vii) {sinkx siiiky)ao. From 
Fig. ^a.) the interorbital operators (v) and (vi) have a 
much larger GS overlap than the intraorbital operator 

(vii) . The mean- field calculation of the gaps for the sym- 
metric combination of the prevailing B2g pairing opera- 
tors, i.e. (v)-|-(vi), is in Fig.|31[c). All the gaps have nodes 
along the x and y axes, also in good agreement with [23| . 
A strong k dependence is observed for all FS pockets, 
and the electron-pocket gaps arc small 1 meV). 

Dynamical Pair Susceptibilities. To complete our 
analysis the dynamical pair susceptibilities, defined by 
P(uj) = dte"^\Ai,/{t)Al^^{0)), were also studied in 
the state with A^=16 for the pairing operators Ak,j;. A 
procedure used in the context of the cuprates will be 
followed [H. Results for U=3 \ti\, Jh/C/=0.2, and sev- 
eral values of Jnn/U were obtained along the dot-dashed 
line (red) of Fig. [Ifb). The overlaps calculation already 
indicated that for iV=18 there are several low-lying en- 
ergy states with different symmetries near the GS. The 
dynamical pair susceptibilities show that most of these 
low lying states have a large overlap with A]^ j|4'Ar=i6(0)) 

for A]^ with the appropriate symmetry. This is quali- 
tatively different to the cuprates' t-J model, where the 
overlap of the doped GS with Aj^ il^'(O)) was large for A 
with d-wave symmetry but negligible for s-wave symme- 
try In that s-wave case the spectral weight in P{uj) 
accumulates at high energies, while P{uj) for the d-wave 
pairing operator showed a well defined sharp peak at the 
GS energy of the doped state This is not the case 
for the two-orbital model. For example, in Fig. El^d) at 
>/NN/f^=0.095, where the doped GS has symmetry i?2g, 
a sharp peak occurs in P{uj) for the i?2g pairing oper- 
ator (v), but a similar behavior is found in P{uj) for 
the Aig pairing operator (i) (the low-lying peak origi- 
nates in a low-lying excited state with Aig symmetry). 
In addition, the susceptibility for a pairing operator (viii) 
(cosfcj; -I- cosfcy)cr3, NN version of the Big operator (ix) 
(cosfcx cosfcy)cr3, is also competitive [2^ (Fig. [3ljd)). 

Conclusions. The effects of NN and NNN Heisenberg 
terms on the symmetry and the binding energy of two 
electrons added to the undoped state of the two-orbital 
Hubbard model were studied using Lanczos techniques 
on small clusters. Quasi-nodal Aig bound states are 
stabilized for physical values of Ja/U, in the interme- 
diate/large U region, in agreement with RPA results 
Our results also indicate that a competing B2g state may 
become stable in physically relevant regimes of C//|<i|. 
In addition, the pairing susceptibility presents low-lying 



excitations with B2g, ^ig, and Big symmetries. Thus, 
pairing correlations with any of these symmetries could 
be stabilized by small modifications in the model param- 
eters, in agreement with 0-0 This suggests that 
a similar sensitivity to small details may occur among 
different compounds of the pnictide family. 
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